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ABSTRACT: 

Simultaneous measurements of ocean microstructure and sound phase shift 
from a stable platform in BASS STRAIT, AUSTRALIA, have provided new relations 
between the statistics of the medium and the statistics of the local sound 
phase speed near the sea surface in the open ocean. Because of dispersion due 
to ambient bubbles, average phase speeds in the frequency range 15 to 100 kHz 
differ as much as 2.5m/sec from the accepted 3MHz "precision" velocimeter 
values down to depths of 7.6m in the presence of wind speeds of 25-30 knots. 
These differential speeds imply average bubble volume fractions of the order 
of 10-7 w ith standard deviations approximately one-fifth of the mean value. 

The differential sound speed is now shown to increase approximately propor- 
tional to the wind speed. The third power decrease of differential speed with 
increasing depth is roughly verified. Under these experimental conditions 
the predominant cause of the local phase fluctuations at 24.4 and 95.6 kHz is 
shown to be bubble activity rather than temperature fluctuations. At 24.4 kHz 
the activity is the random change of number of bubbles. At a frequency such 
as 95.6 kHz, where there is a large resonant bubble population, the predomi- 
nant part of the frequency spectrum of the sound phase modulation is shown to 
be caused by changing bubble radius due to the fluctuating ocean surface wave 
height. The sound phase spectrum mimics the wind wave spectrum given by 
Pierson and Moskovitz to two octaves beyond the frequency of the peak energy, 
at which point the surface pressure effect has dropped low enough for temper- 
ature fluctuations to take over. A theory is presented for prediction of 
these micros tructural sound phase fluctuations from a knowledge of the surface 
wave height spectrum. 
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I. INTRODUCTION 



The objective of this research was to study the relations between 
acoustic phase shifts and the fluctuating parameters of the near sur- 
face ocean. Time series data of the acoustic phase were used to cal- 
culate the dependence of the standard deviation and the mean value of 
the sound speed on frequency, range and depth. Simultaneous measurements 
were made of the fluctuating temperature, two components of particle 
velocity and standard acoustic velocimeter output. 

The research was conducted using the BHP/ESSO Ocean Drilling Plat- 
form "Kingfish A M , Lat. 38°, 39.5* S, Long. 148°, 08. 7 1 E, Bass Strait, 
southeast of Victoria, Australia, as a test platform. The platform is 
about 50 miles from shore and stands in water approximately 250 feet 
deep. Data were collected during two visits to the tower: BASS I 

(20-25 November 1972) and BASS II (12-18 February 1973). 

This report is being issued jointly with the Royal Australian Naval 
Research Laboratory, Garden Island, 2000, Australia. 

II. BACKGROUND 

Following the early work of Minnaert (1933) both Germany'*’ and the 

2 

United States and its allies conducted extensive theoretical and applied 

research during and following World War II to evaluate the importance of 

ocean bubbles in the propagation of underwater sound. By the end of the 

1950s, the physical understanding of bubble action in the presence of a 

3 

sound wave was secure and in the 1960s research spread to many adjacent 

areas including chemical engineering, nuclear physics, and oceanography 

4 

and even meteorology . 
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The question of whether there are bubbles in the ocean, without the 
intrusion of man-made vehicles, has remained controversial. Those who tend 
to doubt the existence of ambient bubbles in the sea quote the laboratory 
experiments of Liebermann^, who demonstrated the rapid loss of bubbles due 
to buoyancy and diffusion. On the other side of the controversy, there have 
been measurements^ showing that bubbles are produced by wind waves and decay- 
ing matter on the ocean floor^. Marine biologists have no doubts about the 
continual generation of bubbles in the sea by photosynthesis, and marine 
life. The root question, of course, is whether those bubbles that are pro- 
duced by physical and biological activities in the sea are so numerous that 
in - situ measurements will always detect ambient bubbles and, more directly 
for the underwater acoustician, the question is whether ambient bubbles are 
numerous enough to influence sound propagation. 

The Ocean Physics Laboratory at the Naval Postgraduate School conducted 

a student research program to determine the excess attenuation and scatter 

8 9 

from near surface bubbles from 1963 through 1965. This early work y was 
summarized in a recent publication"^ which showed the results of photograph- 
ing bubbles and of in - situ acoustical attenuation and scattering measurements 
at the Naval Underwater Center's oceanographic tower in San Diego, California, 
and in Monterey Bay. Measurements were made from the surface to depths of 
45 ft. Bubble volume fractions of the order of 10 or greater in a one 

micron radius increment were typical. The dominant sound frequencies that 
would be affected were in the measurement range, above about 25 kHz, but theory 
indicates that significant dispersion of sound speed should also be found for 
frequencies lower than the measured ones. 

The Ocean Physics Laboratory set out to measure, directly, the change of 
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sound speed with sound frequency, implied by the previous determinations of 
attenuation and scatter near the sea surface. Results of measurements^ from 
the NUC tower at San Diego in October 1971, showed unmistakable evidence of 
dispersion of the speed of sound. Deviations of the propagation speed with 
frequency were of the order of meters/second greater or less than the "theo- 
retical" speed calculated from the traditional parameters of temperature, 
salinity and depth. The standard deviations of the speed measurements were 
much less than the deviations of the measured speed from the theoretical 
speed. Bubble populations were identifiable at depths 4.3m and 9.3m, which, 
when converted to surface values on the basis of constant mass and isothermal 
transport, represent 62 micron (52 kHz) bubbles and 100 micron (32 kHz) 
bubbles. The power spectral density of the time-varying sound speed peaked 
at frequencies less than 0.5 Hz, suggesting that bubbles were entrained at 
orbital frequencies associated with the ocean surface wave system. 

During 1972, our group mounted three simultaneous attacks in order to 

understand better the dynamics of bubble populations in the upper ocean. Trans- 
12 

port analysis of the effects of buoyancy, gas diffusion and vertical velo- 
cities of the medium showed that bubbles which originate at the surface would 
require large downward convective velocities to persist at depth; large 
bubbles which originate in the sediment rise to the surface with little evi- 
dence of gas diffusion; an exponentially decaying volume source (maximum at 
the surface) would show ambient distributions that mimic the source function, 
if the bubbles are large, and would show much greater variation from the 

source distribution as a result of gas diffusion, if the bubbles are small. 

13 

A study of the statistics of fluctuations of sound speed due to bubble 
populations showed that the maximum fluctuations due to changes in number of 
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bubbles would occur at frequencies displaced from the resonant frequency by 

amounts f / ( 1 ± 6) 2 . The standard deviation would be greater at the lower 
R 

frequency than at the upper one. f R = resonant frequency, 6 = damping constant. 

14 

An experimental study of sound speed at the Naval Underwater Center 

oceanographic tower in June 1972, covered the extended frequency range 15 to 

150 kHz in increments of approximately 3.5 kHz at a depth of 3.1 ± 0.3 meters. 

Both average speed and standard deviation of speed were determined at each 

13 

frequency. An analysis of the experimental data showed that the statisti- 
cal theory predictions of frequencies of maximum fluctuation of sound speed 
agreed with the experimental measurements for at least four identifiable 
bubble populations of surface equivalent radii 54, 76, 105 and 200 microns, 
resonant at surface equivalent frequencies of 60, 43, 31, 16 kHz, respec- 
tively. Correlations between sound phase and ocean surface height revealed 
that there are strong fluctuations in the sound speed, in phase with the 
surface height fluctuations, at sound frequencies that correspond to the 
resonant frequencies of the peak bubble concentrations. 

The purpose of the BASS experiments was to repeat, in the waters of the 
South Pacific, some of the San Diego experiments and to look for the statis- 
tics of the correlations with speed fluctuations and its temperature micro- 
structure and medium particle velocity as a function of depth and sound 
frequency. 

III. THEORY 

A. Phase Speed in a Bubbly Medium 
1. Bubbles of One Size Only 

The presence of bubbles in a water medium affects the speed of sound 
because of the changed compressibility. Although there is also a potential 
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effect on the speed caused by the reduced density of the medium, we are 
interested only in the ocean medium, for which the density change can readily 
be shown to be negligible for realistic environments. 

The dependence of the speed on the compressibility and density is given 
by 



The compressibility of the bubbly medium, K, is made up of a part due to the 
water, K Q , and the part involving the change in volume of the bubbles, , 



K - K„ + 

The compressibility of the bubble-free water is most simply expressed in terms 
of the usual value of the speed of sound through water, c Q , and the density, 



p 0 > 



K o = 



^0 °0 



( 2 ) 



The speed through the bubble-free water can be obtained from a simplified form 
of Wilson* s equation 

C 0 = 1449. + 4.6t - 0.055t 2 + 0.017d (2a) 

where t - temperature [°C], d = depth [meters], and salinity is assumed 35 ppt. 
The more sensitive compressibility, due to the bubbles, is given by 




where N = number of bubbles of radius, a 

Av = change in volume for each bubble in incident field, P 

2 ™ 

S = 4TTa - surface area of each bubble 

and = complex magnitude of radial displacement of bubble surface. 
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The compressibility, K, is written in complex notation because of its depen- 

AM*, 

dence on the radial displacement, § , which is itself complex. Simplifica- 

tion occurs if we now insert 5 and also define 

•*\£ 

Z = f R /f - u^/uu (4) 



and 



6 = R/uom 



(5) 



Then, 



K i " 
AA/J 



N4 n a [Z_ ~ 1 - j6] 



P 0 w 2 [(Z 2 



l) 2 * 6 2 ] 



( 6 ) 



and 



K = K + K = 



1 , 4naNlz - 1-16] 

, 2 ,r..2 .X 

0 c o 



P 0 U) 2 C(z - 1) + 6 Z ] 



(7) 



The expression for the speed of sound in the bubbly medium is now 



= (—¥ = 

\p K/ 



\p 






r + j 

|_ w 2 V (Z 2 _ I) 2 + 6 2 



( 8 ) 



and the real part of c is 

r « AAjJ 



Relc} = c 1 — 

.O'O O 



2n a Nc Q 2 (Z 2 - 1) 



00 



(Z 2 - l) 2 + 6 2 



(9) 



Before interpreting our result, it is useful to put the bubble measurement 
in terms of the fraction of gas in bubble form, U(a), defined for bubbles 
of radius, a, by 

U (a) = [N(a)][jTTa 3 ] 

Then, 



Re{c} = c 



1 - 



3UZ 



z 2 -i 



0 2 . 2 / 2 .\2 , .2 

2a k R '^Z — 1J + o 



( 10 ) 
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It is of interest to look at the special cases for frequencies, (a) much 
less than the resonant frequency of the bubbles, (b) much greater than the 
resonant frequency, and (c) close to the resonant frequency. For simplic- 
ity, the reduction is carried out for standard, sea level conditions. 



(a) f « f R , Z » 1 



c = 




(ID 



then 



-2 



, since at sea level ak R = 1.36 x 10 , the speed reduces to 



c = C o [l “ 8.1 x 10 U(a)] and the value is found to be less than c Q , inde- 

pendent of the particular frequency, f, and dependent only on the fractional 
volume of gas to water, U(a). 



(b) f » f R , Z « 1 

At sea level, the speed has the value 

c = c |l + 8.1 X 10 3 U(a) — - — tA (12) 

°L i + 6 2 j 

which is greater than the bubble-free value c Q , dependent on the particular 
frequency and damping constant, reaching the value c Q only as the frequency 
approaches infinity. 



(c) f « f R , Z « 1 

The speed in the bubbly water is c Q at the resonant frequency. 

It rises rapidly to a peak at frequencies slightly greater than f^, and drops 
rapidly to a minimum at frequencies slightly less than f . 

The detailed behavior of the speed, near the resonant frequency, is 
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shown in Fig. 1, which has been drawn for two cases of different numbers of 
bubbles. The change in speed has been presented in terms of the shift in 
speed 6c = (c — c 0 ) at frequency ratio f/f R . It is interesting to observe 
that the maxima and minima of the resonant swings occur at frequencies given 
approximately by f = f R (l ± 6 R /2). That is, for 6 R = 0.10, the peak and 
trough of the resonance are 6 R /2 ~ 57 0 above and below the resonant frequency. 
The magnitudes of the speed at these points are above and below the speed 
in bubble-free water by an amount approximately times the value at the 

zero frequency limit given by equation (11). 

Study of equation (10), which led us to Figure 1, shows that the speed 
of sound is significantly altered even at frequencies far from the reson- 
ance frequency. For example, for the case of 6 R = 0.1, the speed at 

3 

f/f = 2 is still approximately 0.33 (8.1 x 10 Uc ) above the bubble-free 
R o 

value which exists as f approaches infinity. And at f/f = 0.5, the lowered 

K 

3 

value of the speed is 0.33 (8.1 x 10 Uc q ) less than the asymptotic value 
achieved as the frequency approaches zero. 



2. Bubbles of many Sizes. 

The generalization of our results to a cluster of bubbles of 
various radii is accomplished by including additional terms in equation (7). 
It is easy to handle the situation for a broad range of radii, such as 
exists in the ocean, by replacing N(a) by n(a)da = number of bubbles per 
unit volume with radius between a and a + da. Similarly, U(a) is replaced 
by u(a)da to represent the volume fraction of gas in bubble form to liquid. 
Then, instead of K , as in equation (7), we write 



Jb- 




4tt a.n. (a.)da. 
i i v i 7 i 

p uu 2 




(13) 
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The dispersion of sound speed in a bubbly region of water. The abscissa is the relative 
frequency, f/f^, where f^ is the bubble resonance frequency. The ordinate is the relative 
change in sound speed, compared to the bubble-free value. The lower curve shows the stan- 
dard deviation of the phase speed as a result of change of bubble resonant frequency, 
a~^(f R ), and change of number of bubbles, 0“ : (n(a)). 



Because all contributions to the compressibility are very small quantities, 
they add linearly and the speed of sound in the bubble region can be written 



Re{c) 




( 14 ) 



The effect of a mixture of bubble sizes around a peak population at reson- 
ance frequency, f , is to smear the dispersion curve for the speed so that, 
although the magnitudes of the deviation from the bubble free value are 
increased, the frequency range between the peak and the trough is also in- 
creased beyond the value that we calculated in Section A for the single 
bubble size. 



B. Statistics of Phase Speed in a Fluctuating Medium 

The near-surface region of the ocean is a medium that shows spatial and 
temporal fluctuations of the temperature, salinity, pressure, and three com- 
ponents of particle velocity. Sound propagation takes place in such a 
medium only at the cost of sound phase and amplitude fluctuations along the 

acoustic path. The phenomena have been studied theoretically under rather 

16 

major simplifying assumptions . Recent experimental research has extended 
these studies to the near-surface region of the ocean^, in order to search 
for the real world correlations between the sound fluctuations and the 
fluctuations of the medium. 

Fluctuations in the speed of sound, measured^ near the ocean surface, 
have been shown to have a Gaussian probability density function except in 
the regions of large differential speed of sound, where the PDF was double- 
peaked. The temporal correlation time for phase was about the same as for 
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the oceanographic variables. The power spectral density of the time-varying 
phase showed evidence of bubbles entrained at orbital frequencies associated 
with the surface wave system. 

Work at the Naval Postgraduate School‘d has also revealed that orbital 
motion due to surface wave action shifts inhomogeneities of temperature, 
salinity, and bubbles, and causes phase changes of sound propagated near the 
surface. Phase modulation envelopes were found to have smaller correlation 
times and spatial correlation lengths in the mixed layer than in the thermo- 
cline. Correlation times were greater for propagation in the vertical than in 
the horizontal directions in the thermocline. Anisotropy in the thermocline 
was also indicated by the variance of phase fluctuation being greater for 
horizontal than vertical propagation. 

14 

Cross-correlations between ocean surface height and sound phase have 
been shown to be very large at the sound frequencies where dispersion of the 
speed indicated large populations of bubbles. 

For the purpose of this report it will not be necessary to look at all 
of the details of the mechanisms for sound speed fluctuations. The problem 
can be discussed adequately if we simply point out that the variance of the 
speed fluctuation can, in principle, be due to any of the medium descriptors 
mentioned in the first paragraph of this section. However, realistically, 
fluctuations in salinity and density can readily be shown (from eq (2a) to 
be negligible. Furthermore, it is possible to subtract out the effects of 
motion of the medium either by a double path, as in the case of commercial 
velocimeters , or by subtracting out the effect of the motion indicated by a 
velocity flowmeter, as we have done. It is the latter technique that has 
been used in this experiment. This leaves fluctuating bubbles and temperature 
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as the major sources of sound speed variation. 

It is possible to calculate the variance of the speed fluctuations due 
to temperature, simply and effectively enough for the purposes of this study, 
by linearizing equation (2a) and rewriting it as 



c 0 c ref 



(15) 



For water of temperature 19° C, the linearization parameters are: 
c c = 1471.25 m/sec, 

re l = 2.66 (“cr 1 



The variable effect of bubbles is due to (a) changes of total gas in 

bubble form, U, or (b) changes of the bubble distribution with radius, or 

(c) perturbation of the resonant frequency of a predominant population at 

the frequency of the study. The source of fluctuation (a) has been con- 
13 

sidered carefully and it has been shown that the maxima occur at frequencies 
given by f = f^(l ± 6^/2) for resonant frequency, f^, with damping constant, 

V 

There is little information about preferential sizes of bubbles for 
specific bubble-making processes. However, there is no reason to assume that 
the average distribution of bubble numbers with radius at a given depth and 
time of day will change, except by rather small amounts during a propagation 
experiment. We will assume that fluctuation source (b) ,is negligible. 

Finally, (c) perturbation of the resonant frequency can have a profound 
effect on the variation of the speed of sound and it appears that the change 
of water pressure (due to surface wave action) at the level of the experiment 
can modulate the resonant frequency as given by the equation 



f 



R 



1 

2na 




(16) 
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where 



P = ambient pressure 
o 

P q ° ambient density 

and thereby cause the speed of sound to move up and down the steep part of 
the speed curve, as shown in Figure 1. If this action occurs during propa- 
gation of a frequency that corresponds to a major bubble population, a very 
large variance of sound phase can occur. 

The effect of bubbles in changing the speed of sound can also be lin- 
earized if we restrict our interest to one frequency at a time. Then, in a 
bubbly medium 

c = c 0 [l + 6c (U , f ) ] (17) 

where 6c(U,f) is the increment in speed due to bubbles. Combining equa- 
tions (15) and (17), since both bT and 6c(U,f) are very small compared to 
unity, the expression for the speed of sound can be written 

c c r + bT 4- 6c (18) 

ref v 7 

from which the variance of the speed is 

Var {c} = b^ Var i t} + Var {6c} (19) 

Experimental values of the variance of 6c can be used in two ways: 

(a) At very low frequencies, equation (11) leads to 

9c 2 

Var {6c} = — 9 p y Var {u} , f « f R , (20) 

so that the variance of changes of total fraction of bubble volume, Var {u}, 
can be calculated. 

(b) At frequencies that show 1007 o cross-correlation of sound 
phase with the surface wave height, it can be assumed that the high 
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correlation is completely due to resonant frequency fluctuation caused by 
pressure and radius perturbations defined by equation (16) near Z = f^/f = 1- 
A change of f R at constant f can be understood by considering the speed 
change shown in Figure 1 . At f/f R = 1, there is a maximum rate of change 
of c with change of the ratio. For bubbles beyond the peak and trough of 
the dispersion curve, the effect of change of f/f is very much less, and 
goes in the opposite direction. A positive 1007 o cross-correlation with 
wave height would suggest that phase changes near Z = 1 override changes in 
the opposite direction in the outer regions of the dispersion curves of the 
much large and much smaller bubbles. 

We consider the region near Z = 1 by defining 

Z = fjf = 1 + A where A « 1 (21) 

K 



A = 



V 3yp o / p 
2 tt a f 



1 for Z « 1 



( 22 ) 



For the slow changes accompanying the ocean wave system we can assume the 
isothermal relation 



da/dP 



a 

3P ' 



Differentiation of the above equation then gives the rate of change 



dA 

dP 



_5_ 

6P 



(1 + A) 



(23) 



Integration between limits 



*i - z i 



1=0 at Pl - P s (l + ^s) 



and 



*2 



= A 



at P 0 = P 
2 s 



(i + 4iS + 



10m 



h,m \ 

10m/ 
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yields 




h,m/10m 
d t m 

10mJ 



1 + 



where P g = surface ambient pressure - 1 atmosphere 
d,m = experiment depth in meters 
h,m - wave height in meters 



( 24 ) 



The expansion for small increments is used 



1 4- x ) = x — ~ x^ + - - - 



where x 



h,m/10m 

’ 1 + laH 

10m 



to give 



A 

6 



[ h ,m/ 10m 1 

i + ia 
10m J 



(25) 



Using the approximation Z = 1 + A and taking the variance of both sides of 
equation (10), we get 



Var {6c (A)} = Var i c. 



1 - 



Upc Q 2 A 



YP 0 (4A 2 + 6 2 



to order A « 1. (26) 



2 2 

If we assume that 4A « 6 (e.g., much less than a 57 0 variation in f / f 

K 



for a 60 micron bubble) and use equation (24), we get 



Var {6c(A)3 = Var < 



0 



1 - 



0.83 



h,m 

10m 



, 2 2.2 .. , ckm 

6 < 1 + ToS 



(27) 



or 



°-6c <A> ‘ 



0.83 Uc 






h/10 



( 28 ) 



where = rms value of fluctuations in sound speed due to 

perturbations of bubble resonant frequency. 
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o - the dimensionless fraction, rms height of surface wave 

in meters/10 meters. 

Since we are dealing here with very slowly varying surface heights 
(F < 1. Hz) the change of resonant frequency and radius of the bubble will 
occur with negligible phase lag. Therefore, equation (25) can be viewed 
from the point of view of a sinusoidal change of height, h,m = H,m sin fit 
where fi = 2 ttF 

and F = a component of the surface wave spectrum causing a 
simultaneous change of A in equation (25) and of the differential sound speed 
without phase change at the same driving frequency. In the case of a spec- 
trum of surface height, equation (28) can be rewritten 



6c(fi)dfi = 



0.83 Uc 



(v) 2 ( 



n H(Q)dn 

~ , d ,m \c 2 

1 + 



(28a) 



C . Theoretical Basis for the Experimental Measurements 

The physical form of the acoustic sound path is shown in Figures 3 and 
4. Sound from the piston-like transducer propagates through the frame and 
the sound pressure is measured by a series of very small identical 11 hydro- 
phones along the path. For a plane wave, or in the far field of the piston, 
the phase shift along the distance, x, between two hydrophones is measured 
by counting the integral number of wave lengths, n, and determining the 
residual phase angle, 0, which represents the instantaneous fractional wave 
length, 0/360. Then the instantaneous speed is calculated by the product of 
frequency by wave length, using 

C = n + 0/360 (29) 

The assumptions in the derivation, as well as the accuracy of the component 
measurements, are discussed in reference 11. 
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IV. EXPERIMENTAL DETAILS 



A photograph of Kingfish A Tower is presented in Figure 2. A rectan- 
gular aluminum pipe frame of height six feet and length eighteen feet, was 
used to mount the transducers, hydrophones, flowmeter, velocimeter and ther- 
mister (Figs. 3 and 4). To prevent displacement of the hydrophones, the 
hydrophone mounting wires were put under tension and isolated from the frame 
by rubber shock mounts. To minimize any reflections from the frame, selected 
segments were covered with acoustic absorbent rubber (SOAB). See also ref- 
erences 11 and 14. 

The rectangular frame was hung so that the sound source and receivers 
were in a horizontal line. The method of attachment to a nominal 45° steel 
cable is shown in Figure 4. The steel cable had one end held to the ocean 
floor by an anchor, the other end was fixed to the tower. The steel cable 
and polypropylene attachment lines can be seen in Figure 2 and to the right 
in Figure 3 as they enter the water. 

At the center of one end of the frame, pointing in the horizontal (X) 
direction, an USNRD F-33 piston transducer was mounted. On the acoustic 
axis of this transducer, four Atlantic Research LC-10 hydrophones were 
placed at source separations shown in Table 1. 

TABLE 1. 

SOURCE TO HYDROPHONE SEPARATIONS (METERS) 

BASS I BASS II 



LC 


10-1 


1.047 


1.065 


LC 


10-2 


0.998 


0.988 


LC 


10-3 


2.002 


1.990 


LC 


10-4 


1.355 


0.986 
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Fig 2. Photograph of the platform, Kingfish A, from which the experiments 
were conducted. 
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Fig 4. Sketch of equipment frame and components. 





The transmitting instrumentation (Figure 5) consisted of a General Radio 
coherent decade frequency synthesizer, type 1162-A, generating a stable sinus 
oidal output of two volts rms, selectable in frequency. The signal was ampli 
fied by a Hewlett-Packard Power amplifier, 467A, to voltage levels consistent 
with undistorted signals through the receiving instrumentation, and then 
impressed across the transducer. 

The receiving instrumentation (Figure 5) consisted of four LC-10 hydro- 
phones with flat (-llOdB ± 3 dB re 1 volt /microbar) receiving responses from 
10 kHz to 150 kHz. The signals from the hydrophones were amplified 30 dB 
by NUS pre-amplifiers (model 2010-30) before passing through 180 feet of 
water-proof shielded cable. Two selected signals were then bandpass filtered 
by digitally tuned variable Krohn-Hite model 3322 filters (down 24 dB/octave 
outside of band) set with the center frequency equal to the transmitted 
frequency ± 300 Hz. The filtered signals were then fed into a Dranetz model 
305 phase meter. The phase meter measured the phase angle between the two 
ac voltages and provided a dc voltage output (10 ms response time) propor- 
tional to the difference of the phase angles at 10 mv per degree. The dc 
voltage was FM recorded at 3.75 ips on a Precision Instrument model 2101, 
seven track magnetic tape recorder. 

The thermister was mounted on the frame midway between the first two 
hydrophones. The dc voltage output was FM recorded at 0.5 v per degree C. 

The reference sound speed was obtained by an NUS model 2030 - 103 - 12 
velocimeter operating at 3MHz transmit frequency. The frequency output was 
time averaged and manually recorded. The frequency output was also demodu- 
lated to produce a time varying dc voltage which was FM recorded at 20.48 
mv per m/sec sound speed. 
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Fig 5. Block diagram of the electronics for the measurement 
of sound phase, temperature, sound velocimeter output, 
and two horizontal components of material flow velocity. 



Two orthogonal horizontal components of the water particle velocity 
were determined by the mounting of an EPCO water current meter such that the 
velocities measured corresponded to the direction of acoustic propagation 
of the USNRD F-33 transducer and the transverse or Y direction (BASS II). 

These two outputs were also FM recorded at 1.640 volts per M/sec. 

V. DATA PROCESSING 

Data processing was essentially the same as described in previous Naval 
Postgraduate School reports (Refs. 11,14). See Figure 5. 

A series of ten minute and four minute individual records makes up the 
analog recorded time varying sound phase difference, water current components, 
temperature, and velocimeter data. The approach used was to digitize the 
analog data and use a computer to carry out the statistical correlation and 
spectral analysis . 

The recorded data were digitized using the Royal Australian Naval Research 
Laboratory* s PDP-11 computer. The sampling frequency was 6.4 Hz, giving a 
Nyquist frequency of 3.2 Hz. The sampling interval was 0.15625 seconds and 
the frequency resolution was 0.05 Hz. The analog signals were low-passed 
filtered at 3.2 Hz, scaled and written on magnetic disk. The horizontal beam 
(X) direction component of water current was scaled to meters per second. 

For each instantaneous value of phase, the instantaneous value for the speed 
of sound was calculated, in meters per second, and instantaneous value for 
water current in the X direction was algebraically subtracted to correct for 
apparent sound velocity due to the motion of the medium with respect to the 
sound field. 

The mean values, variances and standard deviations were then determined. 
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The mean speed was then compared with the mean speed values given by 
the NUS Velocimeter. 

The temporal autocorrelations were computed and plotted using a maximum 
time lag of 20 seconds. The temporal cross correlation was also computed 
between sound speed and water cross current, and sound speed and temperature. 
From these, the auto-spectrum was computed. 

VI « ANALYSIS OF RESULTS 
A. Average Speed of Sound as a Function of Depth 

The equipment used for the phase comparison determination of the speed 
of sound is capable of yielding values with an absolute error of less than 
one meter per second‘d, provided that corrections are made for ocean currents 
and the small phase difference between the hydrophone-preamplifier electrical 
responses. Column 3 of Table 2 presents the values of the average sound 
speed (c) determined by our phase comparison method, corrected for ocean cur- 
rent, at two frequencies, 24.4 and 95.6 kHz (BASS II). Hydrophone separa- 
tion was 0.988 m. Column 4 presents the concurrently read and averaged 
values of the mean speed (c 0 ) given by the sound velocimeter. In column 5, 
we take the difference between column 3 and 4 in order to arrive at the 
differential velocity, Ac - (c) — (c 0 ). These values are due to bubbles 
plus a correction for the previously mentioned hydrophone electrical phase 
difference. A complete frequency run at depths great enough to validate the 
assumption of a no-bubble region at 13.7 m was not concluded because of the 
failure of the tape recorder. However, from the evidence of the standard 
deviation of sound speeds shown in Table 3, it is quite convincing that the 
depth of 13.7 m can be used to represent the bubble-free medium, and. 
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therefore, to correct for the hydrophone-preamplifier electrical phase shift. 
This has been done to arrive at column 6 of Table 2, the corrected differen- 
tial speed of sound due to bubbles, 6 c = Ac^ — Ac^ 7 * Finally, column 7 
gives the average bubble fraction, U, calculated from equation (11), using 
average values. 

The frequency 24.4 kHz was selected because it was the lowest frequency 
available with the 35 dB S/N ratio required for accurate phase information. 
The values of the corrected differential speed become more negative (down to 
6.1 M) as depth is increased, presumably because there is a peak bubble popu- 
lation just above 24 kHz at surface level which increases its resonant fre- 
quency to higher values with increasing depth. A strong population sweeping 
by the sound frequency of 24.4 kHz would cause larger negative values of 6 c. 
Finally, however, by 7.6 m depth, the decreasing numbers of bubbles becomes 
the dominant factor and the speed approaches the no-bubble value. 

The frequency 95.6 kHz presents a continuously decreasing value of 6 c 
as depth increases. Reference 14 shows that there is a trough of sound speed 
at this frequency, caused by the larger number of bubbles at higher resonant 
frequency. As the depth is increased, existing large bubbles have increas- 
ingly high resonant frequencies which produce a positive increment to add to 
the speed of sound, thereby reducing the negative corrected fifferential 
speed as the depth is increased. The smaller number of bubbles at greater 
depths acts, also, to bring the value of 6 c closer to zero. 

The ratio of the average change in sound speed ascribable to bubbles, 

6 c, to the average speed in non-bubbly water (c^) (equivalent to the veloci- 
meter reading at 3 megahertz) is plotted in Figure 6 as a function of the 
pressure parameter, (3 = 1.33(1 + 0.1 z). In reference 15 it was shown that 
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P 



133 

1 + 0.12 



FIG 6 Variation of the relative change of speed due to 
bubbles, 6c/ (c ) as a function of the pressure parameter, 

P = 1.33/ (1+0. i Z) where Z - depth in meters. The data of this 
experiment, at 24.4 and 95.6 kHz, measured during winds speeds 
of 25-30 knots, are compared with previous attenuation derived 
data obtained during wind speed 6 knots. 
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the ratio 6c/(c 0 ) behaved as (3 . That deduction was based on attenuation 

measurements performed at wind speeds of 6 knots. The slope is plotted in 
Figure 6 for comparison with the present data. It appears that both the 
24.4 and 95.6 kHz data approximately follow that same depth dependence 
except for the contrary slope due to the previously mentioned 24.4 kHz 
resonant bubbles nearer the surface. It is interesting, also, to compare 
the magnitudes of the bubble effect, as a function of wind speed. This is 
possible if one assumes that the waters in BASS Strait are biologically 
similar to those off San Diego, and that the cause of the sound speed dif- 
ference is bubble presence due to the greatly different wind speed. The ratio 
of the differential sound speeds is approximately 6 to 1, which is close to 
the wind speed ratio of the two experiments (28 knots to 6 knots). Therefore, 
the proportionality of volume fraction and differential sound speed to wind 
speed, assumed in reference 15, is roughly verified here. More data of this 
type are needed to determine whether the simple proportionality is indeed 
correct and, if so, the limit of this proposed law. 

B. Average Speed of Sound as a Function of Sound Frequency 

The dispersive effect of ambient bubbles at sea has been predicted on 

the basis of excess attenuation measurements^, and has been measured directly 

11 14 

in two previous experiments * . Although an attempt to repeat the ambient 

experiment in BASS Strait was frustrated by the tape recorder failure in 
BASS II, data of special interest with respect to bubbles near structures was 
obtained during BASS I. The previously discussed electronic phase shift, 
measured in BASS II, and smoothed between 24 and 96 kHz was used to reduce 
the data. A photograph of the frothy surface of the water next to the tower 
docking area and about 5 meters above the measuring equipment is shown in Fig. 
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Fig 7. Photograph of the ocean surface at the base of the tower during 
the data collection for the graph of Fig 8. The equipment is at 
depth 15 feet, next to the tower. 
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The sound speed dispersion is shown by solid line in Figure 8. The curve of 
reference 14 is shown by dashed line. The peaks and troughs of differential 
speed are the greatest that we have ever determined, 

C. Average Bubble Fraction 

The estimated asymptotic differential speed of -15 meters/sec, of 
Figure 8, when used in conjunction with equation (11) of section 3, leads to 
the value of the average total fractional volume of gas in bubble form to 
water, U - 1.2 x 10 ^ for the region near the tower. This value is approxi- 
mately 30 times the fraction for ambient bubbles found at 3.3 m depth in a 
6 knot wind condition, 10 times the values at depths to 7.6 m in the 28 knot 

wind (Table 2) away from the tower in BASS II, and 6 times the value measured 

4 

directly by microscope near breaking waves at shore . One might suspect that 
the optical measurement, since it is effective only for the larger bubbles, 
Underestimates the total fraction. It is interesting to observe that the 
peaks of sound speed are at frequencies close to those found in previous work 
although the experimental frame and particular hydrophones and preamplifiers 
were different in the different experiments, partly to eliminate the possi- 
bility of instrumental artifacts. 

D, Fluctuations of Bubble Fraction and Sound Speed 

Fluctuations in sound speed may be caused by fluctuations in temperature, 
salinity, depth, fraction of gas in bubble form, resonant frequencies of 
bubbles or distribution of bubbles with radius. From a consideration of 
Wilson's equation for the speed of sound in a bubble-free region, it is clear 
that salinity and depth variations (due to change of sea surface displacement 
with time) are generally quite unimportant compared to the temperature varia- 
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tions. It is of interest, however, to compare the fluctuations of sound 
speed due to temperature variation with those attributable to the different 
bubble effects. This can be done by means of equation (19). 

Table 3 displays the results of calculations of fluctuations during the 
BASS II depth experiment. The experiment was conducted from 1110 to 1255 
(24.4 kHz) and 1255 to 1442 (95.6 kHz) on Saturday, 17 February 1973. The 
equipment was lowered to five depths, 3.0, 4.6, 6.1, 7.6, 13.7 m, at a loca- 
tion about 100 feet from the tower, which was evidently free of tower effects. 
The hydrophone separation was 0.99 m. The wind speed of 25-30 knots measured 
at height 37 m above the ocean surface leads to a predicted rms height of 
waves, CJJ = 1.1 m, assuming the Pierson-Moskovitz spectrum for the fully 

developed sea. This wind had persisted for about 24 hours prior to the 

acoustics experiment. A swell wave of estimated peak to trough height of 2 
meters appeared to be present. Under these conditions the rms variation of 
the average speeds of sound at 24.4 kHz and 95.6 kHz at five depths are given 
in column two of Table 3. Column four presents the calculated rms values of 
c q which is calculated from the variance of the temperature and the constant, 
b, of the linearized equation (15). The variance of the excess speed due to 
bubbles is obtained by subtraction (equation (19)) and its square root, the 
rms value, is presented in column five, "EXPT". It is seen to be greater 
than the rms value due to temperature, (see footnote Table 3) except for 95.6 

kHz at the greatest depth of the experiment, 13.7 m, where we can assume that 

the bubble population is negligible. Since variances add, it is appropriate 

to compare the variance of sound speed due to temperature fluctuations, 

2 

Var c Q = G , as a percentage of the total variance of sound speed fluctua- 
c o 

tions, as shown in column six. The relative variance due to temperature 
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See following page for footnotes and explanation of terms 
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fluctuations is particularly small at the frequency 24.4 kHz, and particu- 
larly large at both frequencies at the greatest depth due to the deficiency 
of bubbles there. The distinction between the two frequencies will be dis- 
cussed later after we have looked at the temporal correlations. 

Finally, the standard deviation of fluctuations of the bubble fraction 
is shown in the last column for the 24.4 kHz sound, where a calculation is 
possibly by using equation (20). The standard deviations, cr~^ , range from 
one-third to one-eighth of the average bubble fractions themselves, U, given 
in Table 2. 

If one has the situation: (a) the variance of sound speed due to temper- 
ature fluctuations is negligible and (b) the variance of sound speed is 
principally due to perturbation of the resonant frequency of bubbles, it is 
possible to use equation (28) to estimate the rms value of the excess speed 
of sound fluctuations. Assumption (a) is verified by column six of Table 3, 
except, of course, for the greatest depth, 13.7 meters. Assumption (b) 
would be applicable only if we are considering a sound frequency at which 
there is a very large population of bubbles. Such a situation can be iden- 
tified either by a large excess absorption, a large excess speed of sound, 
or a large cross correlation between the sound phase and the surface wave 
height, such as was observed in reference 14. The cross correlation test 
is the best criterion to use but, unfortunately, surface wave heights were 
not recorded in this experiment. In section F we will show that at 95.6 kHz 
the phase fluctuation spectrum is highly responsive to the surface wave 
height spectrum. This is pretty good evidence of a strong bubble population 
resonant to 95.6 kHz, so we have proceeded on that assumption. 

The parenthetical inserts of Column 5 of Table 3, presents the results 
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of our attempt to use the simple linear theory (equation (28)) to predict 
a^ c , the rms value of the speed fluctuations due to bubbles for the fre- 
quency 95.6 kHz. Comparison with the experimental values of column 5 shows 
that our predictions are an order of magnitude too large. It is not dif- 
ficult to find the reason for this overestimation. The flaw is that we 
have used our simple theory (equations 21-28) which was developed for a 
single radius population with A < 0.05 in attempting to predict frequency 
swings that were as great as A = 0.16 in a multiradius population in the 
high seas of BASS II. The theory overestimates the magnitude of the speed 

fluctuations because it uses the slope near f - f for the calculation, 

R 

when, in fact, this slope is the maximum slope and is applicable only for 
bubbles of precisely f = f n . For bubbles slightly off resonance and for 
larger excursions of pressure, the changes are very much less. Further- 
more, for very large excursions in pressure, the phase shift becomes 
Counterproductive 11 because the sections of the dispersion curve over the 

hump start to produce phase shifts that cancel the effect near f = f. 

R 

We believe, therefore, that the order-of -magnitude disagreement between the 
experimental and theoretical parts of column 5 are not a contradiction of 
our theory of the probable cause of the speed fluctuations at 95.6 kHz. 
Additional confirmation follows in section F. 

E . Temporal Correlations 

The tentative conclusions about the acoustical effect of bubbles, which 
we obtained by consideration of the variances, are reinforced when we con- 
sider the temporal correlations. Consider Figure 9, which shows the tem- 
poral correlations of the temperature at the four depths, 3.0, 4.6, 6.1 
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AUTOCORRELATION OF TEMPERATURE 
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Fig 9. Temporal autocorrelation of the temperature fluctuations at depths 
3.0, 4.6, 6.1, and 13.7 in. Measurements made during 24.4 kHz sound 
phase study, BASS II. 



and 13.7 m during the 24.4 kHz run. The gross appearance is of an expo- 
nential or Gaussian correlation function with a periodicity of about 6 
seconds (see reference 18) . This corresponds to a narrow band noise 
centered around a peak at 0.16 Hz, which is approximately the frequency at 
which peak energy density is found for the wind wave system that existed 
at the time of the experiment. It is noted that the correlation period is 
a little shorter (higher frequency) for the temperature fluctuations at 
3.0 m depth compared to the greater depths; this would be accounted for by 
filtering action due to depth. That is, at greater depths, the frequency 
spectrum of the particle velocities associated with the surface wave system 
will be shifted toward a lower peak frequency because of the greater attenu- 
ation and, therefore, the more effective removal of high frequencies. Super- 
imposed on the large period variation is a correlation periodicity of about 
1 second. We believe that this fine structure is an artifact produced by 
oscillation of the cable and/or the frame which held all of the equipment. 

The fine structure is repeated in all the records of temperature, components 
of velocity, and phase at 13.7 m. We believe that this fluctuation exists 
also in the phase at other depths but is overwhelmed by the bubble effects 
that we were studying at nearer surface positions. 

Turning now to the temporal autocorrelations of the sound phase, we 
consider Figures 10 and 11, for 24.4 and 95.6 kHz, respectively. In Figure 
10, we see that, at the greatest depth, 13.7 m, the temporal correlation of 
the sound phase at 24.4 kHz is almost a replica of the temperature correla- 
tions of Figure 9. The fine structure appears with the same periodicity of 
about 1 second, and the 6 second periodicity due to motion of the surface 
wave system is clearly seen as well. Recalling our analysis of the 
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Fig 10. Temporal autocorrelation of phase fluctuations for 24.4 kHz 
sound at depths 3.0, 4.6, 6.1, 7.6 and 13.7 m. BASS II. 



variances, this is a case in which the temperature-caused variance of the 
speed is greater than the bubble-caused variance. It is, therefore, not 
surprising to find the sound phase correlation resembling the temperature 
correlation at 24.4 kHz as the temperature microstructure moves with respect 
to the frame. There is only a hint of the fine structure at the depth 7.6 
meters where Table 3 has shown us that the variance of sound speed fluctua- 
tions due to temperature is only 0.2% of the variance due to bubble effects. 
At shallower depths there is even less fine structure because, even though 
the frame moved, the phase fluctuations were due to bubbles rather than to 
temperature blobs. It is our belief that the decorrelation of the sound 
phase at all depths except 13.7 m is principally a measure of the random- 
ness of the strong bubble populations at these nearer-surface positions. 
There is virtually no response to the ocean wave system at the near-surface 
positions in the case of thq 24.4 kHz sound. 

Figure 11, the temporal autocorrelation of sound phase fluctuations at 
frequency 95.6 kHz presents a completely different story. Only at 13.7 m 
is there a hint of the temperature-motion fine structure. The predominant 
general feature is the narrow band noise character of the correlation func- 
tion with a periodicity of approximately 6 seconds, that is, the frequency 
0.16 Hz again, the center of the surface wave energy spectrum. We propose 

that the explanation of these correlation functions, which are so close to 

23 

what is to be expected of a surface wave correlation function , is due to 
the responsiveness of the resonant frequency of predominant bubble popula- 
tions to surface wave pressure variations as explained in Section III, B, 
and as used in Section VI, D, to calculate the theoretical variance in 
Column 5 of Table 3. We do not have the directly measured bubble distribu- 
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Fig 11. Temporal autocorrelations of phase fluctuations for 95.6 kHz 
sound at de ths 3.0, 4.6, 6.1, 7.6 and 13.7 m. BASS II. 



tions to prove absolutely that there was a large bubble population resonant 
at 95,6 kHz at all depths. However, it should be recalled that the resonant 
frequency is proportional to the 5/6 power of the ambient pressure. This 
means that a large bubble population resonant at frequency f Q at depth 3.0 m 
will be resonant at 2f Q at 18 m. Therefdre, it should not be surprising to 
find effects of large bubble populations, resonant at 95.6 kHz, at all depths 
of this experiment conducted in a sea churned up by 25-30 knot winds. In 
fact, there is ample ancillary evident of a large population of bubbles at 
95.6 kHz, as will be seen in Section F. 

Two temporal cross correlations have been plotted in an attempt to 
further characterize the effect of the medium on the sound phase. Figure 12 
shows the cross correlation of flow velocity fluctuations (the axial compo- 
nent of the velocity of the water with respect to the frame and sound field 
axis) and 24.4 kHz sound phase variations. It is observed that the long 
term periodicity of 6 seconds is again found for both 3.0 m and 13.6 m 
depths, but the fine scale fluctuation of approximately 1 second periodicity 
is clearly identifiable only at the greater depth. The fact that the fine 
structure is almost non-existent on the correlogram for 3.0 m is presumably 
due to the very much larger peak frequency orbital flow velocity components 
which will tend to swamp the other temperature-caused fluctuations of sound 
phase variation after normalizing. Figure 13 is a correlogram of the 
temporal correlation of flow with temperature. It is almost identical to 
the previous figure for the same reasons. 

It is observed that the correlation envelopes at 95.6 kHz require much 
more time to decay than their counterparts at 24.4 kHz, particularly the 
shallow depth correlograms at 24.4 kHz. The interpretation that we make is 
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Fig 12. Temporal crosscorrelation of horizontal flow and 24.4 sound 
phase at depths 3.0 and 13.7 m. BASS II. 
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Fig 13. Temporal cross correlation of horizontal flow and temper- 
ature fluctuations during 24.4 kHz study at depths 3.0 and 
13.7 m. BASS II. 



that there is greater variability in the instantaneous total number of 
bubbles relative to their average number than there is variability in 
temperature. Tables 2 and 3 show that CJy/U ranges from about one-third 
to one- tenth as depth increases. Notice that the near-surface 24.4 kHz 
correlograms are down to 1/e of their value in about 1 second, whereas 
the envelope of the temperature-dominated 24.4 kHz correlogram at 13.7 m 
depth is still above 1/e after as much as 7 seconds' delay. If this 
reasoning is correct, the spatial correlation function of number of 
bubbles (if this could be measured) would have a smaller correlation length 
than the spatial correlation function of excess temperature. 

Using the same reasoning as above, we must conclude that the number of 
resonant bubbles at 95.6 kHz, which determines the character of those cor- 
relograms, must be reasonably constant over a 7 second lag time (the maxi- 
mum calculated). Since the 95.6 kHz bubbles are the smallest ones measured 

in this experiment, it appears that the small bubbles retain their identity, 
Aa 

— ^ « .05, and their gross locations, for a time of at least 7 seconds. 

F . Spectra 

^ ♦ 

Data processing of the spectra of 0, T and v have been done with dif- 
ferent time lags and different resolutions and band width in order to clarify 
the spectral implications already deduced from the correlation functions. 

The previously described "fine structure" that is so prominent in the temper- 
ature correlations (Figure 9) and elsewhere, turns out to be an harmonic 
series with 0.5 Hz fundamental and strong, even, harmonics particularly at 

*To avoid confusion with the symbol, u, used for volume fraction, we will 
use v for the horizontal material velocity in the direction of the sound 
beam. 
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1, 2, and 4 Hz. These components are strongest in the material velocity and 
temperature spectra, and we believe that they represent motion of the equip- 
ment frame. Since this oscillating motion is in the presence of much larger 
D.C. velocities (e.g., v - 0.64 m/sec and - 0.064 m/sec in one case - see 
Appendix 3) the even harmonics are emphasized. Further, since the oscillat- 
ing motion of the frame takes the temperature sensors to positions of dif- 
ferent temperature (due to the vertical temperature gradient), the temper- 
atures show similar spectra. In the cases where temperature effects, rather 
than bubble effects , dominate the fine structure is evident also in the 
sound phase fluctuations. Since the fine structure is an artifact that is 
characteristic of our mounting structure rather than the ocean, we will 
not discuss this aspect of the spectra any further. 

A higher resolution (F* = 0.025 Hz), smaller bandwidth (F max =0.6 Hz), 

presentation of the horizontal velocity component of material velocity at 

19 

five depths is shown in Figure 14. Previous study has shown that the 
spectral components of the material velocity in the ocean are well corre- 
lated with, and representative of, the ocean surface height spectrum, 
attenuated and filtered by the depth dependent factor exp(-Kz), where K is 
the wave number of the spectral component of the orbital velocity. The 
velocity spectra of Figure 14, do, indeed, resemble those of an ocean sur- 
face height spectrum (confirmed also in Figure 12 of the BASS II Trials 
Officers Report^) with peak power spectral density (PSD) at the frequency 
0.175 Hz. 

The behavior of the PSD of $ V (F) at the peak frequency (0.175 Hz) has 

*To avoid confusion, we will use the upper case "F n for modulation 
frequency and the lower case "f M for sound carrier frequency. 
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Fig 14. Power spectral densities of fluctuations in the 
horizontal material velocity at 3.0, 4.6, 6.1, 
7.6 and 13.7 m during the 95.6 kHz sound phase 
study. BASS II. 
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been identified and replotted in Figure 15 as a function of depth in order 
to compare with theory. At this frequency, the surface wave number is 
K = 0.116 m 1. If the spectral component is following the attenuation law 
with depth, 20 log-^Q of exp(-Kz) gives a change of level, with depth of 
1.01 dB/m. A best line of this slope has been drawn through the five data 
points available, and it is seen to quite well represent the behavior of the 
observed material velocity component. We conclude that the peak component 
of the underwater particle velocity spectrum in BASS II was principally due 
to the surface wave height spectrum. The point at which the higher frequency 
components of the velocity spectrum are dominated by turbulent (non-orbital) 
components is not clear. 

The spectra of the temperature fluctuations at five depths during the 
95.6 kHz run are plotted in Figure 16. Since the temperature spectrum in- 
volves the manner in which a spatially distributed temperature microstruc- 
ture is convected past the sensor by the complex motion of the medium with 
respect to the sensor, the spectrum is much more difficult to comprehend 
than the velocity spectrum. Nevertheless, some understanding can be deduced. 
For example, the spectra are due to rates of change of the temperature and 
are, therefore, measured by the material change, 



dT 

dt 



ST 

5t 



+ (V 



V)T 



The local rate of change, S/St, can be assumed to be negligible com- 
pared to the convective change if the Taylor Hypothesis of frozen convected 
scalars is accepted. Since vertical temperature gradients are more likely 
than horizontal ones, the convective change will be maximum in the z direc- 
tion, that is, the term w ST/Sz. Unfortunately, only one electromagnetic 
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AT PEAK FREQUENCY 0.17 Hz 




Fig 15. Power spectral densities of horizontal component of 
material velocity fluctuations at the peak frequency 
of 0.17 Hz, "x", plotted as a function of depth, and 
compared with the theoretical attenuation for orbital 
velocities of surface waves, exp(-Kz) = 1.01 dB/m. 
BASS II. 
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Fig 16. Power spectral densities of fluctuations of 
temperature at depths 3.0, 4.6, 6.1, 7.6 and 
13.7 m. Measurements made during 95.6 kHz 
sound phase study. BASS II. 
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flowmeter was available and it was set to measure the horizontal components 
rather than the vertical. However, our previous discussion of Figure 14 
suggests that the component at 0.17 Hz was at the peak of both the pressure 
and the velocity spectrum, and, therefore, that the z component of the 
orbital velocity, w, would also peak at this frequency. Now consider Fig- 
ure 16, We see that even though the vertical component of the material 
velocity is greater at 3.0 than at 4.6 m, the temperature fluctuation spec- 
trum at the peak frequency is less at the shallower position. This can only 
be due to a vertical temperature microstructure that has a large gradient 
dT/Sz at the lower level. This reversal of the expected decrease in $rj> with 
increasing depth does not occur at any other level, at the peak frequency. 

An inversion at approximately 5 m in the BT at 1010 and a "twitch 11 at about 
4 m in the BT at 1415 (see Appendix 2) support our observation that the 
microstructure at 4.6 m was distinctive. The effect also shows in the CHp 
column of Table 3, during the 95.6 kHz run. 

Another feature of the temperature spectra is that the prominent peaks 
of energy at 0.17 and 0.27 Hz at the shallower depths appear to be shifting 
to higher frequencies as the depth is increased. This effect appears also 
in the velocity spectra, and may be evidence of underwater turbulence, uncor- 
related with the surface height, which is becoming increasingly significant 
as the depth increases. If this is the true explanation, it appears that 
the turbulence spectrum has energy peaks at higher frequencies than those 
of the orbital velocity spectrum. There is some confirmation of this specu- 
lation in the greater than expected value of the velocity at 13.6 m in 
Figure 15. The temperature spectrum appears to imitate these peaks in the 
velocity spectrum and therefore to be behaving somewhat like the passive 
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convected scalar of Batchelor* s theory 

Spectra of the fluctuations of sound phase are shown in Figures 17, 18, 
and 19. Compared to 95,6 kHz, the spectra for sound frequency 24,4 kHz are 
generally flatter and closer to a Gaussian form, as would be expected from 
the near-Gaussian correlation functions (figures 17, 19). The broader spectra 
for 24.4 kHz are undoubtedly a consequence of the strong dependence of the 
sound speed on the bubble fraction, U, which is likely to be Gaussian dis- 
tributed (statistically) and to have a Gaussian spatial correlation function. 

We have already pointed out the broad PDF implied by 0.12 ^ ^ 0.30. There 

U 

is essentially no evidence of the peak surface wave frequency (0.17 Hz) in 
the 24.4 kHz spectra. 

The predominant feature of the phase fluctuations for 95.6 kHz is the 
strong peak at the ocean surface wave peak frequency of 0.17 Hz*. Figures 
18 and 19 show this peak to be common to all depths and to generally decrease 
in magnitude with increasing depth; this provides additional evidence of our 
explanation that the fluctuations at this frequency are predominantly due to 
change of radius caused by the surface wave pressure. Away from this domi- 
nant frequency, the spectra at depths 3.0, 6.1, 7.6 and 13.7 m are almost 
indistinguishable. Beyond F = 0.3 Hz, the 95.6 kHz spectrum at depth 4.6 m 
shows larger PSDs than at 3.0 m (as does the phase spectrum for frequency 
24.4 kHz); this is apparently due to the larger temperature fluctuations at 
4.6 m previously observed in Figure 16. Except for this anomaly at 4.6 m, 

*For these graphs the resolution is 0.05 Hz so that what was formerly 
identified as a peak at 0.175 Hz, now is found at 0.15 rather than 0.20 Hz; 
in all cases we have compromised and arbitrarily assigned the unresolvable 
value 0,17 Hz to this peak frequency so that we will not cause confusion in 
the identification by the reader. 
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Fig 17. 



Square root of power spectral densities of 24.4 kHz 
sound phase fluctuations at depths 3.0 and 4.6 m. 
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Fig 18. 



Square root of power spectral densities of 95.6 kHz 
sound phase fluctuations at depths 4.6, 6.1 and 7.6 m. 
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Fig 19. Square root of power spectral densities of 24.4 kHz 
sound phase fluctuations at 3.0 m and 95.6 kHz sound 
phase fluctuations at 3.0 and 13.7 m. BASS II. 
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for F > 0,17 Hz, the power spectral density of the phase fluctuations for 

95.6 kHz sound follows precisely the law that would be expected for a 

21 22 

surface wave pressure effect ’ for a fully-developed sea. This is shown 
most dramatically in Figure 20, where we have replotted some of the data of 

Figures 18 and 19 on log-log paper in order to show the power law dependence. 

At low frequencies the phase fluctuation spectrum is almost a twin to the 
ocean surface wave spectrum and the fifth power law for F > F m is found at 
all depths for a range of frequencies from just past F m = 0.17 Hz to approxi- 
mately F = 0.5 or 0.7 Hz. The exception, again, is the case for depth 4.6 m 

where the fluctuating temperature effect begins to take over from the pres- 
sure dependence for F > 0.35 Hz (graphs for 4.6 m, 6.1 m, and 13.7 m not 
shown to avoid crowding). The displaced lines of slope F ^ have been re- 
plotted versus depth (not shown). The variation of the spectral level with 
depth has been determined to be approximately exp(-H/4.7 m) , which compares 
fairly well with a determination of (n/a)da from decay of attenuation with 
depth (reference 10), that behaved as exp(-H/6.7 m) . 

As a conclusive check to show that the spectrum of phase fluctuations 
at 95.6 kHz was not determined by the particle velocity, the phase shift due 
to spectrum of underwater material velocity components, $(F) - Ae ^ F 
was plotted separately and compared with the data for the 3.0 m spectrum of 
phase fluctuations (not shown); the velocity theory did not fit the data. 

We conclude that the phase modulation at 95.6 kHz is due to changes of radius 
of resonant bubbles in response to pressures caused by the major frequency 



components of the surface wave height spectrum. 




Fig 20. Square root of power spectral densities of 95.6 kHz 
sound phase fluctuations. BASS II. 
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VIII. APPENDICES 



APPENDIX 1 



BASS I EXPERIMENTS 



I. Data 

1. Fluctuations experiment at separations 1.047, 2.045, 4.046 m 

2. Decorrelation rate: data at 15, 25, 35 ft depth. One separa- 

tion. 62.15 kHz. 

3. Bubble sweep: data at 28 frequencies from 15 to 97 kHz, 15 ft 

dep th . 

BASE II EXPERIMENTS 



I. Data 

A. Acoustical data 

1. Frequency run: four minute runs from 10 kHz; single wavelength 

differences from 50 to 75 kHz, selected frequencies elsewhere. 
Depths: 10, 70 (max), 20, 10 ft. Estimated Duration of Record- 

ing (EDR) : 4 hours; Estimated Duration of Experiment (EDE) : 

6 hours . 

2. Depth run: 10 minute runs at two frequencies, at six depths. 

a. Set frequency at 25 kHz*, vary depth 10, 15, 20, 25, 50, 
and 70 feet. 

b. Set frequency at 95 kHz, vary depth 70, 50, 25, 20, 15, 
and 10 feet; EDR: 2 hours; EDE: 3 hours 

3. Diurnal run: sound speed as function of time 

a. Set depth 10 feet below trough. Consecutive 10 min runs at 
25 kHz and 65 kHz; repeat every hour for 24 hours. 

EDR: 8 hours; EDE 24 hours 



^Frequency 25 kHz is frequency for which S/N is at least 35 dB 
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4. Wind run: sound speed during changing wind conditions 

a. Set depth 10 feet below trough. Do continual recording at 
25 kHz" during rapidly changing wind conditions. 

EDR: 3 hours; EDE: 3 hours. 

5. Range run; sound speed as a function of range; frequencies 25 

and 95 kHz at four ranges (10 min runs) at depth 10 feet, 
then repeat at 70 feet. 

EDR: 3 hours; EDE: 4 hours 

B. Auxiliary data needed during all NPS experiments 

1. Two horizontal components of flow velocity (u and v) ; temper- 
ature, T; velocimeter fluctuations and frequency averages. 

C. Hourly observations needed during NPS experiments 

1. XBT from tower 

2. Meteorological observations: wind direction, speed, sky 

coverage, precipitation (height of anemometer) 

3. Oceanographic observations: sea state and swell height, and 
periodicity and direction, white caps, foam. 

II. Data Processing 

A. Time series graphs to determine validity of recordings 

B. Mean values of all data 

C. Variance wrt mean value (= zero) for all data 

D. Temporal Correlations: corrected c, all data of depth run (6 depths, 
2 frequencies) and six-10 min extracts of wind run v, T: one depth 
run and six-10 min extracts of wind 

E. Auto spectra: corrected c at 10 and 70 feet; at 25 and 95 kHz plus 

6 winds, u,v,T: at 10 and 70 feet plus u T , only for 6 winds 

F. Cross Spectra, Coherences and Phase Angles: for same data as in F . 

NOTES: Speed of sound, c, referred to above, is always speed of sound 

corrected for phase shift due to longitudinal velocity, u. 
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APPENDIX 2 



BATHY THERMOGRAMS 
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APPENDIX 3 



MATERIAL VELOCITIES CALCULATED FROM ELECTROMAGNETIC 
FLOWMETER OUTPUT DURING BASS II - 17 FEB 1973 



u = horizontal speed in direction of sound path 
v = horizontal speed perpendicular to sound path 



Concurrent 



Clock 

Time 


Depth 

m 


sound freq. 
kHz 


<u> 

m/ sec 


% 

m/sec 


<v> 
m/ sec 


Tv 

m/sec 


1110 


3.0 


24.4 


0.23 


0.15 


0.33 


0.18 


1235 


13.7 


24.4 


0.64 


0.064 


— 


- 


1258 


13.7 


95.6 


0.64 


0.071 


— 


— 


1318 


7.6 


95.6 


0.61 


0.08 


- 


- 


1343 


6.1 


95.6 


0.56 


0.09 


- 


- 


1401 


4.6 


95.6 


0.51 


0.10 


- 


- 


1430 


3.0 


95.6 


0.53 


0.11 


— 


— 



NOTE: v component output failed after 1130. 
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